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We formulate and solve the equations governing the transonic behaviour of a general relativis- 
tic black-hole accretion disc with non-zero advection velocity. We demonstrate that a relativistic 
\^ ' Rankine-Hugoniot shock may form leading to the formation of accretion powered outflow. We show 

f^ , that the critical points of transonic discs generally do not coincide with the corresponding sonic 

f^ ■ points. The collection of such sonic points forms an axisymmetric hypersurface, generators of which 

^SJ ' are the acoustic null geodesies, i.e. the phonon trajectories. Such a surface is shown to be identical 

with an acoustic event horizon. The acoustic surface gravity and the corresponding analogue hori- 
zon temperature Tah at the acoustic horizon are then computed in terms of fundamental accretion 
parameters. Physically, the analogue temperature is associated with the thermal phonon radiation 
analogous to the Hawking radiation of the black-hole horizon. Thus, an axisymmetric black-hole ac- 
cretion disc is established as a natural example of the classical analogue gravity model, for which two 
kinds of horizon exist simultaneously. We have shown that for some values of astrophysically relevant 
accretion parameters, the analogue temperature exceeds the corresponding Hawking temperature. 
We point out that acoustic white holes can also be generated for a multi-transonic black-hole accre- 
tion with a shock. Such a white hole, produced at the shock, is always flanked by two acoustic black 
p^^ ' holes generated at the inner and the outer sonic points. Finally, we discuss possible applications of 

^^ , our work to other astrophysical events which may exhibit analogue effects. 
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I. INTRODUCTION 

A. Black-hole analogue 



o 

"*T* ■ In recent years, strong analogies have been established between the physics of acoustic perturbations in an inhomo- 
j^ , geneous dynamical fluid system, and some kinematic features of space-time in general relativity. An effective metric, 
IV ■ referred to as the 'acoustic metric', which describes the geometry of the manifold in which acoustic perturbations 
. !^ ' propagate, can be constructed. This effective geometry can capture the properties of curved space-time in general 
S^ ' relativity. Physical models constructed utilizing such analogies are called 'analogue gravity models' (for details on 
%-i I analogue gravity models, see, e.g. the review articles [llQj and the book Q). 

One of the most significant effects of analogue gravity is the 'classical black- hole analogue'. Classical black- hole 
analogue effects may be observed when acoustic perturbations (sound waves) propagate through a classical, dissipa- 
tionless, inhomogeneous transonic fluid. Any acoustic perturbation, dragged by a supersonically moving fluid, can 
never escape upstream by penetrating the 'sonic surface'. Such a sonic surface is a collection of transonic points in 
space-time, and can act as a 'trapping' surface for outgoing phonons. Hence, the sonic surface is actually an acoustic 
horizon, which resembles a black-hole event horizon in many ways and is generated at the transonic point in the fluid 
flow. The acoustic horizon is essentially a null hyper surface, generators of which are the acoustic null geodesies, i.e. 
the phonons. The acoustic horizon emits acoustic radiation with quasi thermal phonon spectra, which is analogous 
to the actual Hawking radiation. The temperature of the radiation emitted from the acoustic horizon is referred to 
as the analogue Hawking temperature. 

In his pioneering work, Unruh [^ showed that the scalar flcld describing the acoustic perturbations in a transonic 
barotropic irrotational fluid, satisfied a Klein-Gordon type differential equation for the massless scalar field propagating 
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in curved space-time with a metric that closely resembles the Schwarzschild metric near the horizon. The acoustic 
propagation through a transonic fluid forms an analogue event horizon located at the transonic point. Acoustic 
waves with a quasi thermal spectrum will be emitted from the acoustic horizon and the temperature of such acoustic 
radiation is given by |j] 



^AH = l -J^ , (1) 

47rKB ICs or] j Acoustic Horizon 

where kb is Boltzmann's constant, h is the Planck constant devided by 27r, Cs the speed of sound, u± the component 
of the flow velocity normal to the acoustic horizon and d/drj represents the normal derivative. The temperature Tah 
defined by Eq. ^ is the acoustic analogue of the usual Hawking temperature Th : 

" STTKeGAfbh 

and hence Tah is referred to as the analogue Hawking temperature. In eq. ^, Mi,^ is the black- hole mass, G Newton's 
gravitational constant, c the velocity of light in the vacuum. Note that the sound speed in Eq. (Q in Unruh's original 
treatment |J] was assumed constant in space. 

Unruh's work |j] was followed by other important papers [a, M S B ■ ^ more general treatment of the classical 
analogue radiation for Newtonian fluid was discussed by Visser |3| who considered a general barotropic, inviscid fluid. 
The acoustic metric for a point sink was shown to be conformally related to the Painleve-GuUstrand-Lemaitre form of 
the Schwarzschild metric and a more general expression for analogue temperature was obtained, where unlike Unruh's 
original expression |j|, the speed of sound was allowed to depend on space coordinates. In order to determine the 
analogue Hawking temperature of a classical analogue system, one needs to know the location of the acoustic horizon, 
the velocity of the fluid and the speed of sound and their space gradients at the acoustic horizon. 

In the analogue gravity systems discussed above, the fluid flow is non-relativistic in flat Minkowski space, whereas 
the sound wave propagating through the non-relativistic fluid is coupled to a curved pseudo-Riemannian metric. This 
approach has been extended to relativistic fluids ^9] by incorporating the general relativistic dynamics of the fluid 
flow. Since the introduction of viscosity may destroy Lorenz invariance, the acoustic analogue is best studied in a 
vorticity free dissipationless fluid. 

B. Transonic accretion as a blaclc-hole analogue: the motivation. 

The process by which any gravitating, massive, astrophysical object captures its surrounding fluid is called accretion. 
If Cs{r) is the local speed of sound and u{r) is the instantaneous radial velocity of the accreting fluid, moving along a 
space curve parameterized by r, then the local Mach number of the fluid can be deflned as M{r) — u{r)/cs{r). The 
flow will be locally subsonic or supersonic, according to M{r) < 1 or M{r) > 1, i.e. according to u{r) < Cs{r) or 
u{r) > Cs(r). The flow is transonic if at any moment it crosses AI — 1. This happens when a subsonic to supersonic 
or a supersonic to subsonic transition takes place cither continuously or discontinuously. The points where such 
crossing takes place continuously are called sonic points, and the points of discontinuous transition are called shocks 
or discontinuities. 

If the accreting material is assumed to be at rest far from the black hole, the flow must exhibit transonic behaviour 
in order to satisfy the inner boundary conditions imposed by the event horizon. Since the publication of the seminal 
paper by Bondi in 1952 [l3|, the transonic behaviour of accreting fluid onto compact astrophysical objects has been 
extensively studied in the astrophysics community. Similarly, Unruh's paper Hi initiated a substantial number of 
works in the theory of analogue Hawking effects with diverse fields of application 1, 3]. However, except for Moncrief 
[ll| |74| and Anderson [l^l , until recently no attempt was made to bridge the astrophysical black-hole accretion and 
the theory of analogue Hawking radiation, by providing a self-consistent study of analogue Hawking radiation for real 
astrophysical fluid flows, i.e. by establishing the fact that accreting black holes can be considered as a natural example 
of analogue system. Since both the theory of transonic astrophysical accretion and the theory of analogue Hawking 
radiation are based on general relativity, it is almost self-evident that accreting black holes can be considered as a 
natural example of analogue system. 

Motivated by the above mentioned arguments, it has recently been shown |13L Il4 | that a spherically accreting 
astrophysical black-hole system is a unique example of classical analogue gravity model which exhibits both the black- 
hole event horizon and the analog acoustic horizon. Hence, an accreting astrophysical black hole may be considered 
an ideal candidate to study these two different types of horizons theoretically and to compare their properties. 

Analogue effects in an idealized axisymmetric system have recently been investigated |1.5j. There, the axisymmetric 
astrophysical accretion has been modelled by a rotating fluid disc of constant thickness with the distribution of matter 



in the disc independent of the z coordinate along the axis of rotation. In this paper, we study analogue effects in a 
more realistic astrophysical system in which the accretion disc thickness is calculated by solving the relativistic Euler 
equation in the vertical direction. The expression for the local disc height obtained in this way is a function of the 
local fluid velocity, the local speed of sound and the space coordinate r. As will be shown in section IIIFI in our 
realistic disc model, in contrast to the disc with constant thickness, the sonic points do not coincide with the critical 
points, and the overall picture may also differ from that in J15l |. 

In the following sections, we describe how to model a general relativistic, axially symmetric, multi-transonic flow 
of perfect fluid accreting onto astrophysical black holes. We formulate and solve the basic equations in section ^ 
and then explore the transonic behaviour of the flow in section HTll In section Hvl we show that a relativistic standing 
shock wave can form in the accretion disc if the flow is potentially multi-transonic. We discuss the relevant acoustic 
geometry in detail in section|V| In section IVTl we show how to calculate the analogue Hawking temperature. Finally, 
we conclude the paper with section IVlII 

II. TRANSONIC BLACK-HOLE ACCRETION DISC IN GENERAL RELATIVITY 

A. Multi-transonic flow 

For the flow of matter with non-zero angular momentum density, accretion phenomena are studied employing 
axisymmetric configuration. Accreting matter is thrown into circular orbits around the central accretor, leading to 
the formation of accretion discs. The pioneering contribution to study the properties of general relativistic black- 
hole accretion discs may be attributed to two classic papers pii Il7| . For certain values of the intrinsic angular 
momentum density of accreting material, the number of sonic points, unlike in spherical accretion, may exceed 
one, and accretion is called 'multi-transonic'. The study of multi-transonic flow was initiated by Abramowicz and 
Zurek llSl. Subsequently, multi-transonic behaviour in black-hole accretion discs have been expansively studied 
El mUl iillllilllEmiHIHI^llS^I. Typically, the outermost sonic point lies close to the corresponding 
Bondi radius. The innermost sonic point and the middle sonic point exist within and outside the marginally stable 
orbit, respectively, for the general relativistic as well as for the post-Newtonian model of accretion flow. The location of 
the sonic points can be calculated as a function of the specific flow energy £ (Bernoulli's constant), the specific angular 
momentum A and the inflow polytropic index 7. The literature on multi-transonic flow usually deals with low angular 
momentum accretion flow. Sub-Keplerian[75| weakly rotating flows are exhibited in various physical situations, such 
as detached binary systems fed by accretion from the so-called OB stellar winds 31, j2], semi-detached low-mass 
non-magnetic binaries l33| , and super-massive black holes fed by accretion from slowly rotating central stellar clusters 
|34l l35| . Even for a standard Keplerian accretion disc, turbulence may produce such low angular momentum flow 
(see, e.g. J3g and references therein). 

B. The dynamics 

For the most general description of fluid flow in strong gravity, one needs to solve the equations of motion for the 
fluid and the Einstein equations. The problem may be simplified by assuming the accretion to be non-self-gravitating, 
so that the fluid dynamics may be dealt with in a metric without back-reactions. We use the units G = c = Mth = 1, 
so that radial distances and velocities are scaled in units Vg = 2'''' and c, respectively, and all other derived quantities 
are scaled accordingly. 

We use the Boyer-Lindquist coordinates |33 with signature — 1-++, and an azimuthally Lorentz boosted orthonor- 
mal tetrad basis corotating with the accreting fluid. We define A to be the specific angular momentum of the flow 
and neglect any gravo-magneto-viscous non-alignment between A and the black-hole spin angular momentum. 

Let Vf^ be the four velocity of the (perfect) accreting fluid. The energy momentum tensor T'^" is then given by 

T^,^ ^ {<^+p) VtJ.Vu + P9t^u, (3) 

with e and p being the fluid energy density and pressure, respectively. 

In this paper, we study the inviscid accretion of hydrodynamic fluid. Hence, our calculation will be focused on the 
stationary axisymmetric solution of the energy momentum and baryon number conservation equations 

T^^. = 0; {pvn,^ = 0, (4) 

where p is the rest-mass density. Specifying the metric to be stationary and axially symmetric, the two generators 
^'^ = (d/dt)'^ and 0'^ = {d/dcf)^ of the temporal and axial isometry, respectively, are Killing vectors. 



We consider the flow to be 'advective', i.e. to possess considerable radial three-velocity. The above-mentioned 
advective velocity, which we hereafter denote by u and consider it to be confined on the equatorial plane, is essentially 
the three- velocity component perpendicular to the set of hypersurfaces {5]„} defined by u^ = const, where v is the 
magnitude of the 3-velocity. Each S^ is timelike since its normal 77^ ex d^v^ is spacelike and may be normalized as 
77^77^ = 1. 

We then define the specific angular momentum A and the angular velocity f2 as 

A = -^; »=4-- '*^^^'" , (5) 

vt v^ g^^ + \gt^ 

The metric on the equatorial plane is given by [iTJ 

ds^ = g^^dx^'dx" = —dt'^ + — (d(f)- ujdt) + —dr"^ + dz^ , (6) 

A r'^ A 

where IS. — r^ — 2r + a^ , A — r'^ + r^a^ + 2ra^, and uj = 2ar/A, a being the Kerr parameter related to the black-hole 
spin. The normalization condition v^Vf^ — —1, together with the expressins for A and Q in Eq. |(SJ), provides the 
relationship between the advective velocity u and the temporal component of the four velocity 
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Ar^A 



1/2 

(7) 



(1 - u2) {^2 - AXarA + X^r"^ {Aa^ - r'^A)} 

C. Thermodynamics 

In order to solve Eqs. I@J, we need to specify a realistic equation of state. In this work, we concentrate on polytropic 
accretion. However, polytropic accretion is not the only choice to describe the general relativistic transonic black-hole 
accretion. Equations of state other than the adiabatic one, such as the isothermal equation 25* or the two-temperature 
plasma (3^], have also been used to study the black-hole accretion flow. 

We assume the dynamical in-fall time scale to be short compared with any dissipation time scale during the accretion 
process. To describe the fluid, we use a polytropic equation of state of the form 

P - Kp\ (8) 

where the polytropic index 7 equal to the ratio of the two speciflc heats Cp and Cy of the accreting material is 
assumed to be constant throughout the fluid. A more realistic model of the flow would perhaps require a variable 
polytropic index having a functional dependence on the radial distance, i.e. 7 = 7(7"). However, we have performed 
the calculations for a sufflciently large range of 7 and we believe that all astrophysically relevant polytropic indices 
are covered. 

The constant K in Eq. ^ may be related to the specific entropy of the fluid, provided there is no entropy generation 
during the flow. If in addition to © the Clapeyron equation for an ideal gas holds 

P^^^pT, (9) 

^rUp 

where T is the locally measured temperature, ^ the mean molecular weight, mn^rnp the mass of the hydrogen atom, 
then the speciflc entropy, i.e. the entropy per particle, is given by |4S| | 

1 1 

■ log i^ H — - + constant, (10) 



7- 1 7- 

where the constant depends on the chemical composition of the accreting material. Equation (|10|l conflrms that K 
in Eq. ((HJ is a measure of the specific entropy of the accreting matter. 
The specific enthalpy of the accreting matter can now be defined as 

h^^J^±^, (11) 

P 

where the energy density e includes the rest-mass density and the internal energy and may be written as 

e = p+-^. (12) 

7- 1 



The adiabatic speed of sound is defined by 



From Eq. (IT^ we obtain 



Combination of Eq. I|13() and Eq. (jH)) gives 



dp 



(13) 



constant entropy 



dp /7-1-c^ 



de \ 7 — 1 
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Using the above relations, one obtains the expression for the specific enthalpy 

'*= ] 2 ■ 

7 - 1 - c; 



(14) 



(15) 



(16) 



The rest-mass density p, the pressure p, the temperature T of the flow and the energy density e may be expressed in 
terms of the speed of sound Cg as 
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7 - 1 - c^ 



(17) 



p = K -'-1 



7 - 1 - ci 



(18) 
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(19) 



e^ K T-i 



7- 1\ ^^"' 



7-1-c^ 
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7 V7- 1-c^ 



(20) 



D. Disc geometry and conservation equations 

We assume that the disc has a radius-dependent local thickness H, and its central plane coincides with the equatorial 
plane of the black hole. It is a standard practice in accretion disc theory to use the vertically integrated model in 
describing the black-hole accretion discs where the equations of motion apply to the equatorial plane of the black 
hole, assuming the flow to be in hydrostatic equilibrium in the transverse direction. We follow the same procedure 
here. The flow variables are averaged over the disc height, i.e. a quantity y used in our model is vertically integrated 
over the disc height and averaged as y = /g dhy/H. We follow [33 to derive an expression for the disc height H 
in our flow geometry since the relevant equations in |39j | are non-singular on the horizon and can accommodate both 
the axial and a quasi-spherical flow geometry. In the Newtonian framework, the disc height in vertical equilibrium is 
obtained from the z component of the non-relativistic Euler equation where all the terms involving velocities and the 
higher powers of {{z/r) are neglected. In the case of a general relativistic disc, the vertical pressure gradient in the 
comoving frame is compensated by the tidal gravitational field. We then obtain the disc height 



H 



7+1 



(7-l)c? 



{l-{l + cl)}{X^vl-a^[vt-l)} 



(21) 



which, by making use of Eq. Q, may be be expressed in terms of the advective velocity 



The temporal component of the energy momentum tensor conservation equation leads to the constancy along each 
streamline of the flow specific energy E (relativistic analogue of Bernoulli's constant) defined as [IJ 



£ — hvt 



From Q and ^^ it follows 



(7-1) 



7-(l + c2) 



l-u^ 



Ar^/\ 



A^ ~ AXarA + X^r'^{'ia? - r^A) 



(22) 



(23) 



The rest- mass accretion rate M is obtained by integrating the relativistic continuity equation Q. One finds 

1 "' 



VT 



(24) 



Here, we adopt the sign convention that a positive u corresponds to accretion. The entropy accretion rate S is a 
quasi-constant multiple of the mass accretion rate: 



S=[-\ '" 47rA^cl-'-^^ 



(7-1) 



VT^^h-a + c^) 



(^) 



H. 



(25) 



Note that, in the absence of creation or annihilation of matter, the mass accretion rate is a constant of motion, whereas 
the entropy accretion rate is not. As the expression for S contains the quantity K = p/p^, which measures the specific 
entropy of the flow, the entropy rate S remains constant throughout the flow only if the entropy per particle remains 
locally unchanged. This latter condition may be violated if the accretion is accompanied by a shock. Thus, 5 is a 
constant of motion for shock-free polytropic accretion and becomes discontinuous (increases) at the shock location, 
if a shock forms in the accretion. One can solve the two conservation equations for £ and S to obtain the complete 
accretion profile. In this paper, we concentrate on the Schwarzschild metric only. A more general solution for the 
Kerr metric is in progress and will be presented elsewhere. For a — 0, the expressions for £, M and S are 



^s, 



chwarzschild 



(7-1) 
7-(l + c2) 
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(26) 



-^^^Schwaizschild — 



Anpcsr^u /2 (7 - 1) [r^ - A^ (r - 2)] 
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(27) 
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The corresponding disc height is given by 



r4-AV(r-2)]' 
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Schwarzschild ^^ 
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7[7-(l + c2)](r-2) 



(28) 



(29) 



E. Velocity gradients and critical points 



By taking the logarithmic derivative of both sides of Eq. H28() we obtain the sound speed gradient as 



dcs 
dr 



,(7-l)[7-(l + c?)] 



where 
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u dr 



r4~ A2r(r-2) 



(30) 



(31) 



DifFerentiation of both sides of Eq. H26|) and the substitution of {dcg/dr) from Eq. (|30|) gives the advective velocity 
gradient 



du _ (:^)c^/i(r,A)-/2(r,A) _M{r,X,c,) 
dr " ,. 
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where 



/2(r,A) 



2r 
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2r3 - aV + A" 

r4 - X^r{r - 2) 



(32) 



(33) 



A real physical transonic flow must be smooth everywhere, except possibly at a shock. Hence, if the denominator 
V {u, Cs) of Eq. H32|l vanishes at a point, the numerator J\f (r, A, Cs) must also vanish at that point to ensure the physical 
continuity of the flow. One therefore arrives at the critical point conditions by making V {u, Cg) and Af (r, A, Cs) of Eq. 
(|32|l simultaneously equal to zero. We thus obtain the critical point conditions as 



= ±1 



/2(r-c,A) 



7+1 



/i(rc,A) + /2(r„A)^ 



/2(rc,A) 



/i(^c,A) 



(34) 



where Uc = u{rc) and Cc = Cs{rc), re being the location of the critical point or the so-called 'fixed point' of the 
differential equation H32|l . Hereafter, a subscript c marks a quantity evaluated at r = r^. The + or — sign in H34() 
corresponds to accretion or wind, respectively. 

Hereafter, we denote by [£, A, 7] our three-parameter space characterizing the flow behaviour. It is understood that 
[£, A]C[£, A, 7] will be a sub-set of [£, A,7] for a fixed value of 7. A point in the parameter space, with particular 
values of the parameters, is denoted by (f , A, 7) or (£, A). 

The conserved specific energy can be expressed in terms of re 



£ 



2/i(rc,A)(7-l) 



2/i(rc,A)(7-l) + /2(re,A)(7 + l) 



/ rg (r, ~ 2) {/i(r„ A) (7 - 1) + /2(rc, A) (7 - 1)} 
/i(r„A)(7-l){r3-A2(r,-2)} 



(35) 



For a particular value of (£, A, 7), one can now solve Eq. (|35|l to find the corresponding value of r^- 

To determine the behaviour of the solution in the neighbourhood of the critical point, we need to evaluate the 
space gradient of the advective velocity u' = du/dr at the critical point. Equation H32|l is equivalent to the set of two 
parametric first-order differential equations 



-^ =N{r,\,Cs)., 



(36) 



dr 
d7 



V{u,Cs) , 



(37) 



with A/" and T) vanishing simultaneously at the critical or fixed point. The value uj, that u' takes at the fixed point is 
obtained by applying the L' Hospital's rule to Eq. 132f) . We obtain a quadratic equation 



^u'e + 3?i< + 3fi2 = 0, 



with the solutions 



where 
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FIG. 1: Classification scheme for critical points. For details see section Hi El 



The solutions are classified according to their topological behaviour in the neighbourhood of the fixed points. If 5fti 
is non-zero, the inequality 3?i^ < 43^2^ indicates that the fixed point is of a spiral type, and if Ki^ > 43fi2^, the fixed 
point is of a saddle type for 9^2 > or a nodal type for 3^2 < 0. If ^i — 0, there are no nodal points, and the spiral 
points become of a circular type when ^2 < 0. The saddle-type fixed points in this case correspond to the condition 
^2 > 0. The classification scheme for the critical points is depicted by the flow chart diagram shown in Fig. ^ 

For a particular value of {£, A, 7), Eqs. H26|l - (|43|l can be solved simultaneously to obtain a transonic solution passing 
through the critical point represented by (£", A,7). However, one should note that an 'acceptable' physical transonic 
solution must be globally consistent, i.e. it must connect the radial infinity r^oo with the black-hole event horizon 
r = 2rg. This acceptability constraint further demands that the critical point corresponding to the fiow should 
be of a saddle or a nodal type. This condition is necessary although not sufficient, as we discuss in the following 
sections. Spiral or centre-type critical points do not admit a globally acceptable solution. For further discussion on 
the classification of critical points, see, eg. |2a.l40l | and references therein. 



F. Distinction between critical points and sonic points 

Before we proceed, let us discuss one important issue regarding the relation between the critical points and the 
sonic points in this work. From Eq. (|34|) one can calculate the Mach number of the flow at the critical point as 



Clearly, M^ is generally not equal to one, and for 7 > f, is always less than one. Hence we distinguish a sonic point 
from a critical point. In the literature on transonic black-hole accretion discs, the concepts of critical and sonic points 
are often made synonymous by defining an 'effective' sound speed leading to the 'effective' Mach number (for further 
details, see, eg. [19|,|2^). Such definitions were proposed as effects of a specific disc geometry. We, however, prefer to 
maintain the usual definition of the Mach number for two reasons. 

First, in the existing literature on transonic disc accretion, the Mach number at the critical point turns out to be a 
function of 7 only, and hence Mc remains constant if 7 is constant. For example, using the Paczyiiski and Wiita pdl | 
pseudo-Schwarzschild potential to describe the accretion phenomena leads to 



However, the quantity Mc in Eq. (|44(l is clearly a function of r^ and hence, generally, it takes different values for 
different re for multi-transonic accretion, even at a fixed value of (£, A, 7). In the following paragraphs we show that 
the difference between the radii of the critical point and the sonic point may be quite significant. We define the radial 
difference as 

Ar^ = |r, -re|. (46) 

The quantity Ar* may be a complicated function of (£, A,7), the form of which can not be expressed analytically. 
The radius Ts in Eq. 1)46(1 is the radius of the sonic point corresponding to the same (£, A,7) for which the radius 
of the critical point r,. is evaluated. Note, however, that since r^ is calculated by integrating the flow from Tc, Ar^ 
is defined only for saddle-type critical points. This is because, as we will see in the subsequent sections, a physically 
acceptable transonic solution can be constructed only through a saddle- type critical point. One can then show that 
Ar* can be as large as 10^ rg or even more (for details, see section UlI E|l . 

The second and perhaps the more important reason for keeping Tc and Ts distinct is the following. In addition to 
studying the dynamics of general relativistic transonic black-hole accretion, we are also interested in studying the 
analogue Hawking effects for such accretion flow. We need to identify the location of the acoustic horizon as a radial 
distance at which the Mach equals to one, hence, a sonic point, and not a critical point will be of our particular 
interest. To this end, we first calculate the critical point for a particular (£,A,7) following the procedure discussed 
above, and then wc compute the location of the sonic point by integrating the flow equations starting from the critical 
points. The details of this procedure are provided in the following sections. Furthermore, the definition of the acoustic 
metric in terms of the sound speed does not seem to be mathematically consistent with the idea of an 'effective' sound 
speed, irrespective of whether one deals with the Newtonian, post Newtonian, or a relativistic description of the 
accretion disc. Hence, we do not adopt the idea of identifying critical with sonic points. However, for saddle-type 
critical points, r^ and r^ should always have one-to-one correspondence, in the sense that every critical point that 
allows a steady solution to pass through it is accompanied by a sonic point, generally at a different radial distance r. 

It is worth emphasizing that the distinction between critical and sonic points is a direct manifestation of the non- 
trivial functional dependence of the disc thickness on the fluid velocity, the sound speed and the radial distance. In 
the simplest idealized case when the disc thickness is assumed to be constant, one would expect no distinction between 
critical and sonic points. In this case, as has been demonstrated for a thin disc accretion onto the Kerr black hole 
|15| . the quantity Ar^ vanishes identically for any astrophysically relevant value of 5, A, 7, and the Kerr black-hole 
spin parameter a. 

III. GLOBAL CLASSIFICATION OF THE [E,\,'y] SPACE BASED ON THE NATURE OF THE 

CRITICAL POINTS 

A. Choice of [£,7] and the classification scheme 

Next, we give a complete classification scheme for the critical points and the corresponding flow solutions, for the 
parameter space spanned by all astrophysically relevant values of [£,A,7]. We first set the appropriate bounds on 



10 

[£",7] to model the realistic situations encountered in astrophysics. Since the specific energy £ includes the rest-mass 
energy, £" = 1 is the lower bound which corresponds to a flow with zero thermal energy at infinity. Hence, the 
values £" < 1, corresponding to the negative energy accretion states, would be allowed if a mechanism for a radiative 
extraction of the rest-mass energy existed. The possibility of such an extraction would in turn imply viscosity or other 
dissipative mechanisms in the fluid, the properties which would violate Lorenz invariance. Since Lorenz invariance is 
a prerequisite for studying the analogue Hawking effects, we concentrate only on non-dissipative flows, and hence, we 
exclude £ < 1. On the other hand, although almost all £ > 1 are theoretically allowed, large values of £ represent flows 
starting from infinity with very high thermal energy. In particular, £ > 2 accretion represents enormously hot flow 
configurations at very large distance from the black hole, which are not properly conceivable in realistic astrophysical 
situations. Hence, we set 1 ^ f ^ 2. 

The physical lower bound on the poly tropic index is 7 = 1, which corresponds to isothermal accretion where 
accreting fluid remains optically thin. Hence, the values 7 < 1 are not realistic in accretion astrophysics. On the 
other hand, 7 > 2 is possible only for superdense matter with a very large magnetic field and a direction-dependent 
anisotropic pressure. The presence of a magnetic field would in turn require solving the general relativistic magneto- 
hydrodynamic equations, which is beyond the scope of this paper. Thus, we set 1 ;$ 7 ^ 2. However, astrophysically 
preferred values of 7 for realistic black-hole accretion range from 4/3 (ultra-relativistic) to 5/3 (purely non-relativistic 
flow) 42]. Hence, we mainly focus on the parameter range 



4 5' 

1<£:<2, -<7<- 



(47) 



Figure 121 shows a complete classification scheme for the critical point parameter space together with some represen- 
tative flow topologies. The figure has been drawn for the fluid with 7 — 4/3, but a similar figure may be produced for 
any other value of astrophysically relevant 7. The central plot of Fig. |21 depicts the classification of [f , A] according 
to the number and nature of the critical points in the flow. The side panels show the representative topologies for all 
distinct critical point zones drawn in the central plot. The panel marked by, say, W shows a representative topology 
for the particular (£, A) chosen from the region marked by W in the central plot. A thick black cross with a white 
square at its centre in the W zone indicates the point in the parameter space (£, A) for which the panel topology W 
is drawn. The horizontal axis in the panel figures represents the radial distance in units of Vg on the logarithmic scale, 
whereas the vertical axis represents the flow Mach number. The parameters £ and A for which a particular topology 
is drawn are written at the top right corner of the corresponding panel figure. 

B. Mono-transonic solutions with one critical point 

The wedge shaped regions marked by O and I in Fig. |21 represent the values of [£, A] for which there exists only one 
critical point, and hence only one sonic point. The accretion is mono-transonic and the critical points are of a saddle 
type. In the region marked by I, the critical points are called inner critical points since these points are quite close 
to the event horizon, approximately in the range 2 < r""'^'<10. In the region marked by O, the critical points are 
called 'outer critical points', because these points are located considerably far away from the black hole. Depending 
on the value of (£, A, 7), an outer critical point may be as far as lO^rg, or more. The corresponding flow topologies 
are shown in the panel diagram marked by I and O. 

To describe the procedure for obtaining the panel plots, first consider the topology I. Using the value of {£, A, 7) 
marked in the figure, we first solve Eq. H35() to obtain the corresponding critical point re — 4.535 marked in the 
figure by B at the intersection of the accretion branch ABC and the wind branch DBE. We then calculate the critical 
value of the advective velocity gradient at r^ from Eq. H39l) - (|43|l . By integrating H3UI) and 1)32(1 from the critical point 
B, using the fourth-order Runge-Kutta method, we then calculate the local advective velocity, the polytropic sound 
speed, the Mach number, the fluid density, the disc height, the bulk temperature of the flow, and any other relevant 
dynamical and thermodynamic quantity characterizing the flow. In this way we obtain the accretion branch ABC by 
employing the above mentioned procedure. 

Each solution represented in the panel plots in Fig. [3 is two-fold degenerate owing to the ±u degeneracy which 
reflects the physical accretion/wind degeneracy. We have, however, removed the degeneracy by orienting the curves, 
and thus each line represents either the wind or accretion. We have arbitrarily assigned the + sign solution in Eq. 
(|39|l to the accretion and the — sign solution in Eq. H39|l to the 'wind' branch DBE. This wind branch is just a 
mathematical counterpart of the accretion solution (velocity reversal symmetry of accretion), owing to the presence 
of the quadratic term of the dynamical velocity in the equation governing the energy momentum conservation. 

The term 'wind solution' has a historical origin. The solar wind solution first introduced by Parker [43 has the 
same topology profile as that of the wind solution obtained in classical Bondi accretion jlj . Hence the name 'wind 
solution' has been adopted in a more general sense. The wind solution thus represents a hypothetical process, in 
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FIG. 2: Parameter space division for general relativistic accretion and wind. The side panels show the representative topologies 
for all distinct critical point zones drawn in the central plot. The bold capitals that mark the side panels correspond to those 
that mark the regions of the parameter space in the central plot. The points in the parameter space for which the side panels 
are drawn are marked in the central plot by a thick black cross with a white square at its centre. The inward oriented curves 
marked by arrow heads indicate accretion, while those oriented outwards indicate wind. 



which, instead of starting from infinity and heading towards the black hole, the flow generated near the black-hole 
event horizon would fly away from the black hole towards infinity. The topology of such a process is represented by 
the wind solution DBE. 

The above procedure for obtaining the flow topology is also applied to draw the mono-transonic accretion/wind 
branch through the outer sonic point, i.e. the topology marked by O with r""*'''" = 318.63. In fact, the same 
procedure may be used to draw real physical transonic accretion/wind solutions passing through any acceptable 
saddle-type critical point. Note, however, that AB in the topology I or O, is not the complete subsonic branch, 
because B is a critical point and not a sonic point. Using the procedure described above we have to integrate the flow 
from B to the sonic point where the Mach number equals one. The sonic point is basically equivalent to the acoustic 
horizon at which the analogue Hawking radiation is emitted. We discuss this issue in detail in sections and IVII 
Note that the regions I or O in the central plot do not include the 'zero energy' accretion, i.e. a flow parameterized by 
£ — 1, although from the diagram it may so appear. Zero energy mono-transonic accretion does not admit a steady 
solution passing through the inner or the outer critical point. Note also that the parameter space region I extends up 
to f = 2 along the horizontal axis, which is not shown in the figure for convenience. 
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C. Mono-transonic solutions with two critical points 

The region marked by T in the central plot, including the region marked by Ti with [£,A]rj, C [£,A]rp, represents 
[£, A] for which two critical points exist. We call these two critical points rj and r^, with rl < r^. The T region is 
extended along the horizontal axis to the left up to A = 0, and along the vertical axis up to f = 2. A small patch 
of [£,X] It not shown in the figure exists in the range [1.0002<£<1.0014, 5.18<A<5.28]. Except this small patch, 
almost the whole region in the [£, A] space beyond the nearly vertical line passing through A = 4 generally does not 
admit any real solution for r^. The critical points obtained in the region [£, Aj-p can be further classified roughly into 
three regions for £ > 1. 

The first subset [iJ'iAJrp is depicted in the central plot as a dark triangular zone. For this subset, both sets of 
critical points rl and r^ lie within the radial distance of lOr^, with 3.1<r^<3.1014 and 8.6<r^<10. We have found 
that rl is not associated with any steady solution passing through it, whereas there exists a complete mono-transonic 
accretion/wind solution passing through r^ . However, an accretion solution belonging to this class is topologically 
not very different from any mono-transonic solution obtained for (£, A) €[£, A]j or (£, A) g[£, A]q. One such solution is 
shown in the panel figure marked by Ti. The critical point B for such solution is located at Vc = 9.75. By 'complete' 
solution, we refer to a solution that extends from r = 2 to r— >oo. 

The second subset of [£, A]rp corresponds to the cases where r^ ranges from 2 to 10, while r^ is substantially larger. 
Similarly, in this region we do not find any steady solution passing through rl and we find a complete steady, transonic 
solution passing through r^ . Such a solution is shown in the side panel marked by T with the critical point at about 
Tc — 335. 

The third subset of [£, A]rp, as mentioned earlier, is [1.0002<£<1.0014, 5.18<A<5.28]. Here rl ranges approximately 
from 60 to 115 and r^ approximately from 150 to 3500. Again, steady solutions passing through rl do not exist. 
However, the solutions passing through r^ are steady, but not complete, because those solutions form an outbound 
loop round the corresponding rl enclosed in the loop. One such representative solution is shown in the side panel Li. 
The critical point r^ = 486.79 is marked by B and the corresponding rl = 64.08 is marked by a crossed circle (g). 

In the above discussion, we have concentrated on the solutions for which £ > 1. In Fig. 01 we plot re as a function of 
A for the zero energy accretion, i.e. the accretion with £ = 1. The solid line AB represents rj, whereas the dashed line 
represents r^ . We find no steady solutions passing through r^ , whereas there exists a solution passing through rl which 
is smooth but incomplete as it forms an inward bound loop round r^ enclosed in the loop. In the panel figure marked 
by L2, we plot one representative of such topology. The critical point rl =5.1 is marked by B and the corresponding 



r: 
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18.19 is marked by a crossed circle (g). The curve ABCD represents the incomplete accretion (ABC)/wind (DBC) 
solution. Interestingly, this structure is topologically similar to the inward bound loop representing the incomplete 
accretion/wind branch of a multi-transonic accretion which we discuss in the next paragraph. 

Note here that the accretion with two critical points can practically be considered as an example of mono-transonic 
accretion. Since one of the two critical points (either rl or r^) does not belong to any steady flow through it, only the 
other critical point admits a steady flow to pass through it, and hence there will be only one sonic point associated 
with it. 

D. Multi-transonic accretion and wind 

The wedge shaped region marked by A in the central figure corresponds to the multi-transonic accretion. For any 
{£, A) e[£, A]^, the fiow has three critical points The corresponding flow topology is shown in the panel figure marked 
by A. The line ABC1C2D represents a complete transonic solution passing through the saddle-type owier critical point 
B at Tc = 3310.88. The solid circle • marks the outer sonic point at 7-°"*'^'' = 2930.47 where the flow Mach number 
is equal to one. The corresponding wind solution is shown by the line EBF. The crossed circle (g) marks the location 
of the centre-type middle critical point at r™"^'*''' = 18.32, through which a steady solution cannot be constructed. 
The line EIGH represents the incomplete accretion (GIH)/wind (EIG) solution passing through the saddle type inner 
critical point I at re = 5.097. The corresponding inner sonic point r^"™'' (not shown in the flgure) is located at 
r^ = 4.903. Clearly, this loop structure is very similar to the zero energy accretion topology shown in the panel plot 
L2. 

The set [f , A]^^ (or more generally [£,X,j]j^) thus produces doubly degenerate accretion/wind solutions. Such two 

fold degeneracy may be removed by the entropy considerations since the entropy rates S'(rJ.""'"') and S'(r°"*°'') are 
generally not equal. For any {£, A) £[£, A] (A), we find that the entropy rate S evaluated for the complete accretion 
solution passing through the outer critical point is less than that of the rate evaluated for the incomplete accretion/wind 
solution passing through the inner critical point. Since the quantity 5* is a measure of the specific entropy density of 
the flow, the solution passing through r°"*®'' will naturally tend to make a transition to its higher entropy counterpart. 
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FIG. 3: Zero energy {£ = 1) accretion flow with two critical points. The solid line AB represents the set of inner critical points. 
The flow passing through an inner critical point has a closed loop topology as shown in the side panel L2 of Fig. |21 The dashed 
line represents the set of outer critical points. There is no steady solution passing through an outer critical point. 



i.e. the incomplete accretion solution passing through r"™''''. Hence, 
ABC1C2D to increase its entropy accretion rate by an amount 



if there existed a mechanism for the solution 



AS ^ 5(C""') - S{r°''"''), 



(48) 



there would be a transition to the incomplete solution GIH. Such a transition would take place at a radial distance 
somewhere between the radius of the inner sonic point and the radius of the accretion/wind turning point (r = 94.9 
for this case) marked by G. In this way one would obtain a combined accretion solution connecting r—^00 with r — 2 
which includes a part of the accretion solution passing through the inner critical, and hence the inner sonic point. 
One finds that for some specific values of [£, A,7] C[£, A,7]^, a standing Rankine-Hugoniot shock may accomplish 
this task. A supersonic accretion through the outer sonic point r°"*'"' can generate entropy through such a shock 
formation and can join the flow passing through r™"''''. Although two shock locations, indicated by the vertical lines 
through Ci and C2 with a downward arrow, are found, only one of the two shocks is stable. We discuss the details of 
such a shock formation in section Hvl 

The wedge shaped region marked by W in the central figure represents the [£, X] zone for which three critical points, 
the inner, the middle and the outer are also found. However, in contrast to [f,A]j^, the set [£,A]-yy yields solutions 
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for which S{r™'^'^'^) is less than S'(r°"*'"'). Besides, the topological flow profile of these solutions is different. One such 
solution topology is presented in the panel plot marked by W. Here the loop-like structure is formed through the outer 
critical point. Note a similarity with the topology Li. This topology is interpreted in the following way. The flow HB 
passing through the inner critical point B at re = 4.25 is the complete mono-transonic accretion flow, and ABC2C1D is 
its corresponding wind solution. The solution EFG passing through the outer critical point E at 7-°"'°'' = 301.53, with 



the corresponding outer sonic point marked by a solid circle • located at 



— 343.312 represents the incomplete 



accretion (EF)/wind (FG) solution. The initially subsonic wind solution passing through r™"'''' encounters the outer 
sonic point at r-^"™"" = (4.31, and becomes supersonic. However, as S'abc2CiD turns out to be less than S'feGi the 
solution branch ABC2C1D can make a shock transition to join its counter solution FEG and thereby increase the 
entropy accretion rate by the amount AS' = 5'abC2CiD ~ '5'feg- Here, two theoretical shock locations are obtained, 
out of which only one is stable. Hence the set [5, A, 7]-^ corresponds to mono-transonic accretion solutions with 
multi-transonic wind solutions with a shock. 

Besides 7 = 4/3, for which Fig. 13 has been drawn, we could perform a similar classification for any astrophysically 
relevant value of 7 as well. Some characteristic features of [£, A] would be changed as we vary 7. For example, if fmax 
is the maximum value of the energy and if Amax and Amin are the maximum and the minimum values of the angular 
momentum, respectively, for [£, A]^ for a fixed value of 7, then [fmax, Amax, Amin] anti-correlates with 7. Hence, as 
the flow makes a transition from its ultra-relativistic to its purely non-relativistic limit, the area representing [f , A]^ 
decreases. 






FIG. 4: The difference Ar^ between of the sonic points radius Va and the critical point radius r^ as a function of the specific 
energy £ and the specific angular momentum A. The plots represent a mono-transonic flow passing through the inner critical 
point (plot A) and the outer critical point (plot B); a multi-transonic accretion passing through the inner critical point (plot 
C) and the outer critical point (plot D). 
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E. Dependence of Ar° on [£, A, 7] 

Here, we study the dependence of the radial difference between the critical points and the sonic points on fundamental 
accretion parameters. In Fig. 0] the difference Ar^ defined by Eq. (|46|l is plotted as a function of the specific fiow 
energy £ and the specific angular momentum A for a mono-transonic flow passing through the inner critical point (plot 
A) and the outer critical point (plot B), and for a multi-transonic accretion flow passing through the inner critical 
point (plot C) and the outer critical point (plot D). Obviously, the radial difference between the sonic point and the 
critical point may be quite significant, especially for the accretion/wind flow passing through the outer critical and 
sonic points. The quantity Ar^ anti-correlates with A for a general flow, whereas it correlates with £ for a flow passing 
through the inner critical point and anti-correlates with £ for a flow passing through the outer critical point. 

F. Onset of chaotic behaviour in multi-transonic flow 

In the central part of flgure 2 the boundary between the O and the A zone, and between the W and the I zone, 
provides interesting information about the possibility of the onset of chaotic behaviour of the multi-transonic accretion 
flow. According to the theory of dynamical systems (see, e.g. [ij), if a small change of the control parameters leads to 
profound effects on the outcome of the whole system, the system may then be considered to exhibit chaotic behaviour. 
Along the boundary between O and A, an extremely small change of A or f (or 7 if there were a transition across 
the bounding surface between [£, A,7]o and [i£^. A, 7]^) would lead to a sharp transition from a mono-transonic to a 
multi-transonic flow (and vice versa). A similar situation arises for the [f , A].^ < — > [^ , A]j transition. 

The above phenomenon is an example of the bifurcation of critical points due to a slight change in the control 
parameters. It is possible to study the onset of chaotic behaviour even for a stationary system. Bray and Moore [4^, 
for example, studied the spin-glass phase in terms of a T = flxed point, and established a criteria for the onset 
of chaos in such a system by proposing an effective Lyapunov exponent. In a similar way we propose an effective 
Lyapunov exponent for a transonic flow solution. Let ui{r) and U2{r) be the values of the advective velocity at any 
radial distance r for two transonic solutions characterized by very close values of the angular momentum Ai and A2 
respectively , i.e. (Ai — A2) — >0, with £ and 7 kept fixed. Let ui{r + R) and U2{r + R) be the corresponding velocities 
for the above-mentioned transonic flow solutions at the radial distance r + R. We then deflne the effective Lyapunov 
exponent L^fj, so that the following condition holds: 

ui{r + R) - U2{r + R) _ r,,R 



ui{r) - U2{r) 



e^^sfi\ (49) 



The quantity L^ff can also be calculated by varying £ instead of A. Next, one can show that the system will be 
sensitive to the initial boundary conditions if ^e// > 0. The preliminary calculation J4a| indeed indicates such values 
of L^ff > for a transition from mono- to multi-transonic flow, and vice versa. Hence, transonic flow solutions exhibit 
chaotic behaviour. Investigation of chaos in black-hole accretion discs may shed a new light on the explanation of the 
variability mechanism of active galactic nuclei (for details about AGN variability, see, e.g. |43|). However, the details 
of the calculation are beyond the scope of this paper and will be presented elsewhere. 

IV. SHOCK FORMATION IN RELATIVISTIC ACCRETION 
A. A general overview 

Perturbations of various kinds may produce discontinuities in an astrophysical fluid flow. By discontinuity at a surface 
in a fluid flow we understand any discontinuous change of a dynamical or a thermodynamic quantity across the surface. 
The corresponding surface is called a surface of discontinuity. Certain boundary conditions must be satisfled across 
such surfaces and according to these conditions, surfaces of discontinuities are classified into various categories. The 
most important such discontinuities are shock waves or shocks. In an adiabatic fiow of the Newtonian fiuid, the shocks 
obey the following conditions ^ : 



[[H] = o, [[p + pu^]]^Q, 



= 0, (50) 



where [[/]] denotes the discontinuity of / across the surface of discontinuity, i.e. 

[[/]] = /2 - /i, (51) 
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with /2 and /i being the boundary values of the quantity / on the two sides of the surface. Such shock waves are quite 
often generated in various kinds of supersonic astrophysical flows having intrinsic angular momentum, resulting in a 
flow which becomes subsonic. This is because the repulsive centrifugal potential barrier experienced by such flows is 
sufficiently strong to brake the infalling motion and a stationary solution could be introduced only through a shock. 
Rotating, transonic astrophysical fluid flows arc thus believed to be 'prone' to the shock formation phenomena. 

One also expects that a shock formation in black-hole accretion discs might be a general phenomenon because 
shock waves in rotating astrophysical flows potentially provide an important and efflcient mechanism for conversion 
of a significant amount of the gravitational energy into radiation by randomizing the directed infall motion of the 
accreting fiuid. Hence, the shocks play an important role in governing the overall dynamical and radiative processes 
taking place in astrophysical fluids and plasma accreting onto black holes. The study of steady, standing, stationary 
shock waves produced in black hole accretion has acquired a very important status in recent years. For details and 
for an exhaustive list of references see, e.g. |29|. 

Generally, the issue of the formation of steady, standing shock waves in black-hole accretion discs is addressed in two 
different ways. First, one can study the formation of Rankine-Hugoniot shock waves in a polytropic flow. Radiative 
cooling in this type of shock is quite inefficient. No energy is dissipated at the shock and the total specific energy of the 
accreting material is a shock-conserved quantity. Entropy is generated at the shock and the post-shock flow possesses 
a higher entropy accretion rate than its pre-shock counterpart. The flow changes its temperature permanently at the 
shock. Higher post-shock temperature puffs up the post-shock flow and a quasi-spherical, quasi-toroidal centrifugal 
pressure supported region is formed in the inner region of the accretion disc 29J . 

Another class of the shock studies concentrates on the shock formation in isothermal black-hole accretion discs. 
The characteristic features of such shocks are quite different from the non-dissipative shocks discussed above. In 
isothermal shocks, the accretion flow dissipates a part of its energy and entropy at the shock surface to keep the 
post-shock temperature equal to its pre-shock value. This maintains the vertical thickness of the flow exactly the 
same just before and just after the shock is formed. Simultaneous jumps in energy and entropy join the pre-shock 
supersonic fiow to its post-shock subsonic counterpart. For detailed discussion and references see, e.g. |49ll5(l|. 

In this work, the basic equations governing the fiow are the energy and baryon number conservation equations 
which contain no dissipative terms and the fiow is assumed to be inviscid. Hence, the shock which may be produced 
in this way can only be of Rankine-Hugoniot type which conserves energy. The shock thickness must be very small in 
this case, otherwise non-dissipative flows may radiate energy through the upper and the lower boundaries because of 
the presence of strong temperature gradient in between the inner and outer boundaries of the shock thickness. In the 
presence of a shock the flow may experience the following profile. A subsonic fiow starting from infinity first becomes 
supersonic after crossing the outer sonic point and somewhere in between the outer sonic point and the inner sonic 
point the shock transition takes place and forces the solution to jump onto the corresponding subsonic branch. The 
hot and dense post-shock subsonic fiow produced in this way becomes supersonic again after crossing the inner sonic 
point and ultimately dives supersonically into the black hole. A flow heading towards a neutron star can have the 
liberty of undergoing another shock transition after it crosses the inner sonic point [Tg, because the hard surface 
boundary condition of a neutron star by no means prevents the flow from hitting the star surface subsonically. 

B. The relativistic Rankine-Hugoniot conditions 

For the complete general relativistic accretion fiow discussed in this work, the energy momentum tensor T^", the 
four- velocity v^, and the speed of sound Cs may have discontinuities at a hypersurface E with its normal 77^. Using 
the energy momentum conservation and the continuity equation, one has 

[[pv^]]rj^ = 0,[[Tn]r,,^0. (52) 

For a perfect fiuid, one can thus formulate the relativistic Rankine-Hugoniot conditions as 

[[puT,,]] = 0, (53) 

m,if]]^[[{p + e)vtuT.,]]^0, (54) 

[[T,.V''V1] = [[ip + e)u^Tl+p]]=0, (55) 



where Tu = 1/vT— z? is the Lorentz factor. The first two conditions l|53|) and (|54|l are trivially satisfied owing to 
the constancy of the specific energy and mass accretion rate. The constancy of mass accretion yields 







7-(l + c2) 



0. (56) 
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The third Rankine-Hugoniot condition H55|l may now be written as 



1 ^7 — 1\ ^"^ 



-1 j u^j- cl) + , 

7 y V7-l-c2y' I c2(X-u2) 

Simultaneous solution of Eqs. H56|l and 1)57(1 yields the 'shock invariant' quantity 

"^ (7 - cl) + cl 



= 0. 



5^ = 



sU(1-m2) y^ 



(57) 



(58) 



which changes continuously across the shock surface. We also define the shock strength Si and the entropy enhancement 
<d as the ratio of the pre-shock to post-shock Mach numbers {Si = M^/M+), and as the ratio of the post-shock to 
pre-shock entropy accretion rates (8 = S+/S-) of the flow, respectively. 

C. Shock locations in multi-transonic accretion and wind 

The shock location in multi-transonic accretion is found in the following way. Consider the multi-transonic flow 
topology A in Fig. [21 Integrating along BC1C2D, we calculate the shock invariant Sh in addition to u, Cs and M . 
We also calculate Sh while integrating the sector HIG, starting from the inner sonic point up to the point of inflexion 
G. We then determine the radial distance rsh, where the numerical values of Sh, obtained by integrating the two 
different sectors described above, are equal. Generally, for any value of (5, A, 7) allowing shock formation, one finds 
two shock locations marked by Ci (the 'outer' shock - between the outer and the middle sonic points) and C2 (the 
'inner' shock - between the inner and the middle sonic points) in the figure. According to a standard local stability 
analysis [25J, for a multi-transonic accretion one can show that only the shock formed between the middle and the 
outer sonic point is stable. Hence, in the multi-transonic accretion with the topology A, the shock at Ci is stable and 
that at C2 is unstable. Hereafter, whenever we mention the shock location, we refer to the stable shock location only. 

The topology W in Fig. |21 shows the shock formation in the multi-transonic wind. Here the numerical values of Sh 
along the wind solution passing through the inner sonic point (line BC2C1D) are compared with the numerical values 
of Sh along the wind solution passing through the outer sonic point, and the shock locations Ci and C2 for the wind 
are found accordingly. 

D. Geometry of the accretion with shock and generation of accretion-powered outfiow^ 

As a consequence of the shock formation in an accretion, the post-shock flow temperature will also increase abruptly. 
In Fig. El we plot the temperature of the combined accretion flow BCiIH of the topology A presented in Fig. 13 The 
segment AB corresponds to the pre-shock flow (BCi in the topology A of Fig. [Jl, while the segment CD corresponds 
to the post-shock fiow. The vertical segment BC is a discontinuous increase of the flow temperature due to the 
shock formation. The length of BC, which measures the post- to pre-shock temperature ratio T+/T- is, in general, a 
sensitive function of [£", A, 7]. 

In Fig. (HI we present the disc structure obtained by solving Eq. (|29|l for the combined accretion flow BCiIH with 
the topology A in Fig. 13 The point BH represents the black-hole event horizon. The pre- and post-shock regions of 
the disc are clearly distinguished in the figure and show that the post-shock disc puffs up significantly. It is clear from 
Fig. [Slthat the bulk flow temperature will be increased in the post-shock region. Such an increased disc temperature 
may lead to a disc evaporation resulting in the formation of an optically thick halo. Besides, a strong temperature 
enhancement may lead to the formation of thermally driven outflows. The generation of centrifugally driven and 
thermally driven outflows from black-hole accretion discs has been discussed in the post- Newtonian framework |5lLl53 |. 
The post-Newtonian approach may be extended to general relativity using the formalism presented here. 

Owing to the very high radial component of the infall velocity of accreting material close to the black hole, the 
viscous time scale is much larger than the infall time scale. Hence, in the vicinity of the black hole, a rotating inflow 
entering the black hole will have an almost constant specific angular momentum for any moderate viscous stress. This 
angular momentum yields a very strong centrifugal force which increases much faster than the gravitational force. 
These two forces become comparable in size at some specific radial distance. At that point the matter starts piling 
up and produces a boundary layer supported by the centrifugal pressure, which may break the infiow to produce 
the shock. This actually happens not quite at the point where the gravitational and centrifugal forces become equal 
but slightly farther out owing to the thermal pressure. Still closer to the black hole, gravity inevitably wins and 
matter enters the horizon supersonically after passing through a sonic point. The formation of such a layer may be 
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attributed to the shock formation in accreting fluid. The post-shock flow becomes hotter and denser and, for all 
practical purposes, behaves as the stellar atmosphere as far as the formation of outflows is concerned. A part of 
the hot and dense shock-compressed in-flowing material is then 'squirted' as an outflow from the post-shock region. 
Subsonic outflows originating from the puffed up hotter post-shock accretion disc (as shown in the figure) pass through 
the outflow sonic points and reach large distances as in a wind solution. 

The generation of such shock-driven outflows is a reasonable assumption. A calculation describing the change of 
linear momentum of the accreting material in the direction perpendicular to the plane of the disc is beyond the 
scope of the disc model used in this work because the explicit variation of dynamical variables along the Z axis (axis 
perpendicular to the equatorial plane of the disc) cannot be treated analytically. The enormous post-shock thermal 
pressure is capable of providing a substantial amount of 'hard push' to the accreting material against the gravitational 
attraction of the black hole. This 'thermal kick' plays an important role in re-distributing the linear momentum of 
the inflow and generates a non-zero component along the Z direction. In other words, the thermal pressure at the 
post-shock region, being anisotropic in nature, may deflect a part of the inflow perpendicular to the equatorial plane 
of the disc. Recent work shows that |^ such shock-outflow model can be applied to successfully investigate the origin 
and dynamics of the strong X-ray flares emaneting out from our galactic centre. 

E. Dependence of the shock location on [£',A,7] 

We find that the shock location correlates with A. This is obvious because the higher the flow angular momentum, 
the greater the rotational energy content of the flow. As a consequence, the strength of the centrifugal barrier which 
is responsible to break the incoming flow by forming a shock will be higher and the location of such a barrier will be 
farther away from the event horizon. However, the shock location anti-correlates with E and 7. This means that for 
the same £ and A, in the purely non-relativistic flow the shock will form closer to the black hole compared with the 
ultra-relativistic flow. Besides, we find that the shock strength Si anti-correlates with the shock location Tsh, which 
indicates that the closer to the black hole the shock forms , the higher the strength Si and the entropy enhancement 
ratio Q are. The ultra-relativistic flows are supposed to produce the strongest shocks. The reason behind this is also 
easy to understand. The closer to the black hole the shock forms, the higher the available gravitational potential 
energy must be released, and the radial advective velocity required to have a more vigorous shock jump will be larger. 
Besides we note that as the flow gradually approaches its purely non-relativistic limit, the shock may form for lower 
and lower angular momentum, which indicates that for purely non-relativistic accretion, the shock formation may 
take place even for a quasi-spherical flow. However, it is important to mention that a shock formation will be allowed 
not for every (£, A) e[£, A](A) (or (£, A) e[£, A](W)) . Equation lf55|) will be satisfied only for a specific subset of 
[£, A] (A) (or [£, A] (W)), for which a steady, standing shock solution will be found. 

V. ACOUSTIC GEOMETRY, ACOUSTIC HORIZONS, AND PHONON QUANTIZATION 

Before delving into the calculation of Tah for various kinds of sonic horizon generated in the accretion disc, we 
would first like to discuss the relevant features of the acoustic geometry. 

A. Nonrelativistic acoustic geometry 

Let -0 denote the velocity potential describing the fluid flow in Newtonian space-time, i.e. let u = — VV", where u is 
the velocity vector describing the dynamics of a Newtonian fluid. The specific enthalpy ft, of a barotropic Newtonian 
fluid satsfies V/i = (l/p)Vp, where p and p are the density and the pressure of the fluid. One then writes the Euler 
equation as 

-9ti^ + /i+i(VVj)' + $ = 0, (59) 

where $ represents the potential associated with any external driving force. Assuming small fluctuations around some 
steady background p^ ,p^ and tp^ , one can linearize the continuity and the Euler equations and obtain a wave equation 
in the form (for derivation see appendix IX|) 

d^ iF^d,^) = 0, (60) 

where Lp is the first order term in the expansion ip = ip^ + eip + 0{e'^) and the quantity f^"^ in this equation is a 4 x 4 
symmetric matrix defined by (|A14|I . Equation Ht)U|) describes the propagation of the linearized scalar potential ip. The 
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FIG. 5: Bulk temperature profile of a multi-transonic accretion flow for {£,X,^) shown with a shock at rsh ■ The vertical axis 
is the flow temperature in units of 10^" K and the horizontal axis is the radial distance in logarithmic scale. The sector AB 
represents the pre-shock flow, the sector CD the post-shock flow, and the segment BC is the discontinuous increase of the flow 
temperature due to the shock formation. 



function cp represents the low amplitude fluctuations around the steady background (p^jp^, V''^) and thus describes 
the propagation of acoustic perturbances, .i.e. the propagation of sound waves. 

The form of Eq. H()U|) suggests that it may be regarded as a d'Alembert equation in curved space-time geometry. 
In any pseudo-Riemannian manifold the d'Alembertia operator can be expressed as |54| 



D = 



9m/- Is^.l-g^^s. 



where \gfj,u\ is the determinant and g'^'^ is the inverse of the metric g^^,. Next, if one identifies 



one can recast the acoustic wave equation in the form |3 

1 



V— |G^y| 



^HG^G^") d.V - 0, 



(61) 



(62) 



(63) 
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FIG. 6: Disc structure of a multi-transonic accretion flow with a shock at Vsh- The disc height (vertical axis) is plotted as a 
function of the radial distance (horizontal axis). The point BH represents the black-hole event horizon. 



where G^i, is the acoustic metric tensor for the Newtonian fluid. The explicit form of Gp,y is obtained as 



C^fif = p 



-{cl 



(64) 



Thus, the propagation of acoustic perturbation, or the sound wave, embedded in a barotropic, irrotational, non- 
dissipative Newtonian fluid flow may be described by a scalar d'Alembert equation in a curved acoustic geometry. 
The corresponding acoustic metric tensor is a matrix that depends on dynamical and thermodynamic variables 
parameterizing the fluid flow. 

The acoustic metric (I64|) in many aspects resembles a black-hole type geometry in general relativity. For example, 
the notions such as 'ergo region' and 'horizon' may be introduced in full analogy with those of general relativistic 
black holes. For a stationary flow, the time translation Killing vector ^ = d/dt leads to the concept of acoustic ergo 
sphere as a surface at which G^^^^S,'^ changes its sign. The acoustic ergo sphere is the envelop of the acoustic ergo 
region where ^^ is space-like with respect to the acoustic metric. Through the equation G^^S^^S^" = gu = u^ — c^, it is 
obvious that inside the ergo region the fluid is supersonic. The 'acoustic horizon' can be defined as the boundary of 
a region from which acoustic null geodesies or phonons, cannot escape. Alternatively, the acoustic horizon is defined 
as a time like hypersurface defined by the equation 



0, 



(65) 
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where u± is the component of the fluid velocity perpendicular to the acoustic horizon. Hence, any steady supersonic 
flow described in a stationary geometry by a time independent velocity vector field forms an ergo-region, inside which 
the acoustic horizon is generated at those points where the normal component of the fluid velocity is equal to the 
speed of sound. 

B. Acoustic geometry in a curved space-time 

The above formalism may be extended to relativistic fluids in curved space-time background [3| . The propagation 
of acoustic disturbance in a perfect relativistic inviscid irrotational fluid is also described by the wave equation of the 
form (|63|) in which the acoustic metric tensor and its inverse are defined as 



'" " he. 



[g,,u + (1 - cf )wpWi. 



G"" ^ 



he. 



+ (l--)w^w^ 



(66) 



where p and h are, respectively, the rest-mass density and the specific enthalpy of the relativistic fluid, w^ is the 
four-velocity, and g^i^ the background space-time metric. The ergo region is again defined as the region where the 
stationary Killing vector ^ becomes spacelike and the acoustic horizon as a timelike hypersurface the wave velocity of 
which equals the speed of sound at every point. The defining equation for the acoustic horizon is again of the form 
(|65|l in which the three-velocity component perpendicular to the horizon is given by 



u± = 



iv^v^Y 



iv'^Vp.) 



V'^ilfj. 



where rj^ is the unit normal to the horizon. 

It may be shown that, for an axisymmetric flow, the acoustic metric discriminant defined as 
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(67) 



(68) 



vanishes at the acoustic horizon. A supersonic flow is characterized by the condition V > 0, whereas for a subsonic 
flow, V < [13 • According to the classification of Bercelo et al. 56], a transition from a subsonic {V < 0) to a 
supersonic {V > 0) flow is an acoustic blaek hole, whereas a transition from a supersonic to a subsonic flow is an 
acoustic white hole. Hence, for the relativistic disc geometry presented in our paper, one can show that for multi- 
transonic shocked accretion and wind, an acoustic white hole produced at the shock location Vsh is flanked by two 
acoustic black holes produced at the inner and the outer sonic points. The situation is precisely the same as in the 
accretion geometry with a constant thin disc height considered in 15] where a detailed description of the emergence 
of such black or white holes and of the corresponding causal structures is given. 

C. Quantization of phonons and the Hawking effect 

The purpose of this section is to demonstrate how the quantization of phonons in the presence of the acoustic 
horizon yields acoustic Hawking radiation. The acoustic perturbations considered here are classical sound waves or 
phonons that satisfy the massles wave equation in curved background, i.e. the general relativistic analogue of H63|l . 
with the metric G^i/ given by (|65|l . Irrespective of the underlying microscopic structure acoustic perturbations are 
quantized. A precise quantization scheme for an analogue gravity system may be rather involved '5ff|. However, at 
the scales larger than the atomic scales below which a perfect fluid description breaks down, the atomic substructure 
may be neglected and the fleld may be considered elementary. Hence, the quantization proceeds in the same way as 
in the case of a scalar field in curved space [531 with a suitable UV cutoff for the scales below a typical atomic size of 
a few A. 

For our purpose, the most convenient quantization prescription is the Euclidean path integral formulation. Consider 
a 2-|-l-dimensional disc geometry. The equation of motion H63|l with l|66|l follows from the variational principle applied 
to the action functional 



S[ip] 



dtdrd(t>^/^G^"'^f,ip^^ip . 



(69) 



We define the functional integral 
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(70) 
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where Se is the Euclidean action obtained from (|69(l by setting t = ir and continuing the Euchdean time r from 
imaginary to real values. For a field theory at zero temperature, the integral over r extends up to infinity. Here, 
owing to the presence of the acoustic horizon, the integral over r will be cut at the inverse Hawking temperature 2'k/ k 
where k denotes the analogue surface gravity. To illustrate how this happens, consider, for simplicity, a non-rotating 
fluid (u0 = 0) in the Schwarzschild space-time. It may be easily shown that the acoustic metric takes the form 
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where gu = —(1 — 2/r), u — \vr\/y/—gtt, and we have omitted the irrelevant conformal factor p/{hcs)- Using the 
coordinate transformation 
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dr 



we remove the off-diagonal part from 171|) and obtain 
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Next, we evaluate the metric near the acoustic horizon at r = Tg using the expansion in r — Ts 

{r - Ts) 



at first order 
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and making the substitution 
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where R denotes a new radial variable. Neglecting the first term in the square brackets in H73|) and setting t — ir, we 
obtain the Euclidean metric in the form 



dsi = K^R^dT^ 



dR^ + rl 



where 



-gtt 



Hence, the metric near r = rg is the product of the metric on S^ and the Euclidean Rindler space-time 

ds\ ^dR^ +R^d{KTf. 



(76) 



(77) 



(78) 



With the periodic identification t = t + 2tt/ k, the metric (|78|l describes M^ in plane polar coordinates. 

Furthermore, making the substitutions R = e'^^ / k and = y/rg + tt, the Euclidean action takes the form of the 
2-|-l-dimensional free scalar field action at non-zero temperature 



Sf.M = dr dx dy-{d^,Lpf, 

Jo J-oc J-oo ^ 



(79) 



where we have set the upper and lower bounds of the integral over dy to -|-cxd and — cxd, respectively, assuming that Tg 
is sufficiently large. Hence, the functional integral Z in (|70|l is evaluated over the fields •^{x, y, r) that are periodic in 
T with period 27: j k. In this way, the functional Z is just the partition function for a grandcanonical ensemble of free 
bosons at the Hawking temperature Th = k/ {2'kkb)- However, the radiation spectrum will not be exactly thermal 
since we have to cut off the scales below the atomic scale p . The choice of the cutoff and the deviation of the acoustic 
radiation spectrum from the thermal spectrum is cloasely related to the so-called transplanckian problem of Hawking 
radiation J59j. 

In the Newtonian approximation, equation (|77() reduces to the usual non-relativistic expression for the acoustic 
surface gravity (3 
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VI. ANALOGUE TEMPERATURE OF SONIC HORIZONS IN A BLACK-HOLE ACCRETION DISC 



In general, the ergo-sphere and the acoustic horizon do not coincide. However, for some specific stationary geometry 
they do. This is the case, e.g. in the following two examples: 

1. Stationary spherically symmetric configuration where fluid is radially falling into a pointlike drain at the origin. 
Since u = u± everywhere, there will be no distinction between the ergo-sphere and the acoustic horizon. An 
astrophysical example of such a situation is the stationary spherically symmetric Bondi-type accretion |ly| onto 
a Schwarzschild black hole. 

2. Two-dimensional axisymmetric configuration, where the fluid is radially moving towards a drain placed at the 
origin. Since only the radial component of the velocity is non-zero, u = u± everywhere. Hence, for this system, 
the acoustic horizon will coincide with the ergo region. An astrophysical example is an axially symmetric 
accretion with zero angular momentum onto a Schwarzschild black hole or onto a non-rotating neutron star. 

Here we consider the axisymmetric accretion in which the ergo-sphere and the acoustic horizon do not coincide. 
For a transonic accretion disc, the collection of sonic points at a fixed radial distance forms a surface, the generators 
of which are the trajectories of phonons, because at those points u± = Cg. This surface is an acoustic event horizon. 

The analogue surface gravity at the acoustic horizon can be calculated as a function of fundamental accretion 
parameters, i.e. as a function of [£, A, 7]. In complete analogy to the classical black-hole horizon, the acoustic horizon 
emits the analogue Hawking radiation consisting of thermal phonons. The temperature of the analogue Hawking 
radiation Tah is related to the analogue surface gravity as usual 7] . An axisymmetric transonic accretion disc round 
an astrophysical black hole is thus a natural example of analogue gravity model where two types of horizon (the 
gravitational and the acoustic), exist simultaneously. 

For a stationary configuration, the surface gravity can be computed in terms of the Killing vector 

x'^ = e^ + ncj)'' (81) 

that is null at the acoustic horizon. Following the standard procedure [3, |53 one finds that the expression 



^X'^ = iG^^^Vu-^^iGc.pX^xP) 



(82) 



holds at the acoustic horizon. Here the constant k is the surface gravity, the vector 7/^ is the unit normal to the 
acoustic horizon and d/dr] denotes the normal derivative 
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From equation H82II we deduce the magnitude of the surface gravity as [3 
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where r^ denotes the location of the acoustic horizon. 

For the transonic disc geometry described in section ^ we calculate the norm of x^ at the acoustic horizon as 



V^X^l^,^ = r7V('^.-2)(r3-AV,-F2A2). 
The surface gravity is thus obtained as 
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with the corresponding analogue Hawking temperature given by 



Tah — 
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(87) 



where we have used the units in which h — c — 1. 
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It is now easy to calculate Tah for each kind of sonic points (classified according to section IIII|I as a function of 
[£, A, 7]. We now define the quantity r as the ratio of the analogue to the black- hole Hawking temperature 



Ta 



H 



Th 



This ratio turns out to be independent of the black-hole mass. Hence, we are able to compare the properties of the 
analogue and black-hole horizons for an accreting black hole with any mass, ranging from the primordial black holes 
in the early Universe to the supermassive black holes at galactic centres. 

In Fig. [7| we plot the value of r as a function of the specific energy E and the specific angular momentum A of the 
flow. Four representative cases are shown; 

A Mono-transonic flow passing through the single inner type sonic point. The range of (f , A) here corresponds to 
the central region of Fig. |21 marked by I, except that S is extended up to 2. 

B Mono-transonic flow passing through the single outer type sonic point. The range of(£. A) used to obtain the 
result for this region corresponds to the central region of Fig|21 marked by O. 

C Multi-transonic accretion passing through the inner sonic point. The range of (£, A) corresponds to the central 
region of Fig. |21 marked by A 

D Multi-transonic accretion passing through the outer sonic point. The range of (5, A) corresponds to the central 
region of Fig. |21 marked by A 

Similar figures can be drawn for other types of transonic flow for any value of (£, A, 7) allowing a real physical transonic 
solution. 

It is obvious from Fig. |7|that for 7 = 4/3, the temperature Tah asymptotically approaches Th and is never higher 
than Th- However, we observe that r increases non-linearly with 7, and at very high values of 7, e.g. 7 — 1.63, one 
does obtain [£, A] (mainly for high values of £) for which the analogue Hawking temperature exceeds the black-hole 
Hawking temperature. 

In Fig. |S1 we plot the variation of r with A for a fixed £ for two values of 7. Obviously, the ratio of the analogue 
to the black-hole Hawking temperature keeps increasing with 7. 

VII. DISCUSSION 

In this work, we have established that the general relativistic axisymmetric accretion disc round an astrophysical black 
hole can be considered as an example of classical analogue gravity realized in nature. To accomplish this task, we 
have first formulated and solved the equations describing the general relativistic axisymmetric accretion flow round a 
compact object. We then show that such accretion is transonic, and the collection of sonic points forms an acoustic 
horizon. The acoustic horizon is characterized by analogue surface gravity and the corresponding analogue Hawking 
radiation . We have calculated the corresponding analogue temperature as a function of fundamental accretion 
parameters. We have shown that for some region of parameter space spanned by the specific energy £ , the speciflc 
angular momentum A and the adiabatic index of the flow 7, the analogue Hawking temperature may exceed the 
Hawking temperature of the accreting black hole. Generally, this happens for high values of £ and 7. 

Note that the formalism presented in this work is not restricted to analogue effects in the black hole accretion only. 
Hydrodynamic accretion onto any astrophysical object, e.g. a weakly magnetized neutron star, if it happens to be 
transonic, will produce an acoustic horizon, and thus the analogue effects in such systems can be studied using the 
methodology developed here. Thus, we show in general, that a hydrodynamic transonic accretion in astrophysics may 
be regarded as a natural example of analogue gravity model. 

It is important to note that the accreting astrophysical black holes are the only real physical candidates for which 
both the black-hole event hoizon and and the analogue sonic horizon may co-exist. Hence, our application of the 
analogue Hawking effect to the theory of transonic astrophysical accretion may be useful to compare the properties 
of these two types of horizons. 

However, in our work, the analogy has been applied to describe the classical perturbation of the fluid in terms of 
a field satisfying the wave equation in an effective geometry. It is not our aim to provide a formulation by which 
the phonon field generated in this system could be quantized. To accomplish this task, one would need to show that 
the effective action for the acoustic perturbation is equivalent to a field theoretical action in curved space, and the 
corresponding commutation and dispersion relations should directly follow (see, e.g. [Sg). Such considerations are 
beyond the scope of this paper. Also note that for all types of accretion discussed here, the analogue temperature 
Tah is many orders of magnitude lower compared with the fluid temperature of accreting matter. 
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FIG. 7: The ratio of the analogue to Hawking temperature r (vertical axis) as a function of the specific energy of the flow £ 
and the specific angular momentum A. The representative cases shown are: mono-transonic flow passing through the single 
inner type sonic point (A), mono-transonic flow passing through the single outer type sonic point (B), multi-transonic accretion 
passing through the inner sonic point (C), and multi-transonic accretion passing through the outer sonic point (D). 



In constructing a model for the disc height we have followed Abramowicz et al. J39j, although a number of other 
models for the disc height exist in the literature [13, |2a, |23, ISfl IMl . The use of any other disc height model would 
not alter our conclusion that black-hole accretion disc solutions form an important class of analogue gravity models. 
However, the numerical values of Tah and other related quantities would be different for different disc heights. 

In this work, the viscous transport of the angular momentum is not explicitly taken into account. Viscosity, 
however, is quite a subtle issue in studying the analogue effects for disc accretion. Thirty two years after the discovery 
of standard accretion disc theory [13, [IJ , exact modeling of viscous transonic black- hole accretion, including proper 
heating and cooling mechanisms, is still quite an arduous task, even for a Newtonian flow. On the other hand, from 
the analogue model point of view, viscosity is likely to destroy Lorenz invariance, and hence the assumptions behind 
building up an analogue model may not be quite consistent. Nevertheless, extremely large radial velocity close to the 
black hole implies Ti„/ <C Tmsc, where Tm/ and Tmsc are the infall and the viscous time scales, respectively. Large 
radial velocities even at larg er distances are due to the fact that the angular momentum content of the accreting fluid 
is relatively low |63.l64l *65*|. Hence, our assumption of inviscid flow is not unjustified from an astrophysical point of 



view. However, one of the most significant effects of the introduction of viscosity would be the reduction of the angular 
momentum. We found that the location of the sonic points anti-correlates with A, i.e. weakly rotating flow makes the 
dynamical velocity gradient steeper, which indicates that for viscous flow the acoustic horizons will be pushed further 
out and the flow would become supersonic at a larger distance for the same set of other initial boundary conditions. 
In our model, we have performed the computation of Tah for non-rotating black holes. Recently, the spacetime 
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FIG. 8: The ratio of the analogue to Hawking temperature r as a function of specific angular momentum A for for a fixed 
specific energy E and two values of poly tropic index 7. 



geometry on the equatorial slice through a Kerr black hole has been shown to be equivalent to the geometry experienced 
by phonons in a rotating fluid vortex 66] . Since many astrophysical black holes are expected to possess non-zero spin 
(the Kerr parameter a), a clear understanding of the influence of spin on analogue models will be of great importance. 
Our initial calculation indicates that the black-hole spin enhances the analogue effect, i.e. the analogue Hawking 
temperature increases with the Kerr parameter a. 

In connection to the acoustic geometry, one can define an 'anti-trapped surface' to be a hypersurface in which 
the fluid flow will be outward directed with the normal component of the three-velocity greater than the local speed 
of sound. In stationary geometry, an anti-trapped surface will naturally be constructed by the collection of sonic 
points corresponding to a spherically symmetric or axisymmetric transonic wind solution emanating out from an 
astrophysical source. Transonic outflow (wind) is ubiquitous in astrophysics, spanning a wide range from solar/stellar 
winds to large-scale outflows from active galaxies, quasars, galactic micro-quasars and energetic gamma ray bursts 
(GRB). In section Hill we have shown how to identify the critical and the sonic points corresponding to the wind 
solutions. Such a formalism can be useful in studying the transonic properties of outflow from astrophysical sources. 
Hence our formalism presented in this paper can be applied to study the analogue effects in transonic winds as well. 
Recently Kinoshita et al. 69] performed the causality analysis of the spherical GRB outflow using the concept of 
effective acoustic geometry. Such an investigation can be extended into a more robust form by incorporating our 
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work to study the causal structure of the transonic GRB outflows in axisymmetry, i.e. for energetic directed outflow 
originating from a black-hole accretion disc system progenitor. 

In recent years, considerable attention has been focused on the study of gravitational collapse of massive matter 
clump, in particular, on the investigation of the final fate of such collapse (for a review see, e.g. 70]). Goswami and 
Joshi [71^ 1 have studied the role of the equation of state and initial data in determining the final fate of the continual 
spherical collapse of barotropic fiuid in terms of naked singularities and the black-hole formation. It is tempting 
to study the analogue effects in such a collapse model. Since at some stage the velocity of the collapsing fiuid will 
exceed the velocity of local acoustic perturbation one might encounter a sonic horizons at the radial locations of the 
corresponding transonic points in a stationary configuration. One should, however, be careful about the issue that 
many results in analogue models are based on the assumption of a stationary flow, whereas a collapse scenario is a 
full time dependent dynamical process. 

The correspondence between general relativity and analogue gravity has so far been exploited only on a kinematical, 
i.e. geometrical level. The analogue gravity systems lack a proper dynamical scheme, such as Einstein's field equations 
in general relativity and hence the analogy is not complete. A certain progress in this direction has recently been 
made by Cadoni and Mignemi [Z2,l23| who have established a dynamical correspondence between analogue and dilaton 
gravity in 1-1-1 dimensions. We believe that our approach in which an arbitrary background geometry serves as a 
source for fluid dynamics may shed a new light towards a full analogy between general relativity and analogue gravity. 

APPENDIX A: DERIVATION OF THE ACOUSTIC WAVE EQUATION 

This derivation is based in part on references |3,|4g- The continuity and Euler's equations may be expressed as: 

^ + V-(pu) (Al) 



du 



9u ,^ „, ^ 
- + (..V). 



= -Vp + F (A2) 



with F being the sum of all external forces acting on the fluid wich may be expressed in terms of a potential 

F = -pV$, (A3) 

Euler's equation may now be recast in the form 

-^ = M X (V X u) Vp - V ( -u^ -h $ I (A4) 

ot p \2 J 

Next we assume the fluid to be inviscid, irrotational, and barotropic. Introducing the specific enthalpy h such that 

V/i = — (A5) 

P 

and the velocity potential ip for which u — — Vf/', Eq. I|A4|I may be writen as 

-^ + /i+^(VV^)2 + $ = (A6) 

One now linearizes the continuity and Euler's equation around some unperturbed background flow variables po, po, 
ipQ. Introducing 

/O = P0 + epi+O(e^), P^Po + epi+0{e^), 
ip^i^Q + eiJi+Oie'^), h^ho + ehi, (A7) 

from the continuity equation we obtain 

^ + V • (po uo) = 0; ^ + V • (pi uo + Po ui) = 0. (A8) 



28 



Equation ljA5|l implies 



Using this the linearized Euler equation reads 



dh pi 
hi ^Pi-r- = — • 
dp po 



ot 2 at pi 

Re-arrangement of the last equation together with the barotropic assumption yields 

dp dp ,„ , ^ ^ , > 



Substitution of this into the linearized continuity equation gives the sound wave equation 
Next, we define the local speed of sound by 



d_ 
dt 



dp fdijji 



Po Vi/'i - -^ PqVq I — — + vo ■ yipi 



= 0. 



cl = dp/ dp, 
where the partial derivative is taken at constant specific entropy. With help of the 4x4 inatrix 

-/ : -u 

r^p"" 

. -u : {cl-u^) 
where / is the 3x3 identity matrix, one can put Eq. (|A12|) to the form 



(A9) 



(AlO) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 
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